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Abstract 

We explore a possibility of optimization of the method of determination of 
the top quark mass from Mm distribution in semi-leptonic decays t — » blv at 
LHC and a future linear collider (LC). We discover that the systematic and 
statistical errors of Mt determination can be diminished if considering the high 
moments over the distribution. In the case of LHC this allows one to reduce 
more than in twice the errors, and in the case of LC to approach to the precision 
expected at studying the threshold scan of the total cross-section e + e~ — > ti. 

1 Introduction 

The precision determination of the top quark mass is one of the major research 
problems at colliders of next generation [U El El IH IHj • Being a fundamental parameter 
of the Standard Model (SM), the top quark mass is tightly constrained by quantum 
level calculations with other fundamental parameters. This enables one to test the SM 
and/or to select the probable scenario of its extension on the basis of an independent 
M t measurement. 

A considerable progress in this direction is expected at Tevatron and LHC where 
the accuracy of the M t determination is anticipated of about 1-2 GeV lj. At LHC in 
view of the copious production of the top quarks, for the increase of the accuracy the 
decrease of systematic errors is crucial. An analysis of Ref. shows that the method 
most promising from the point of view of optimization of the errors is based on the 
investigation of a distribution over the invariant mass of the observable products of 
semi-leptonic decays t — > bW — > b tv ; more precisely of the isolated lepton C and the 
fi + ji~ pair indicating a J /if) meson produced from the decay of the b quark [B]. In 
this channel one can obtain experimentally very clean final states. Correspondingly, 
the systematic error of the M t measurement can be made low. The evaluation made 
by Monte-Carlo (MC) modelling gives 0.6-0.8 GeV at the statistical error of about 
1 GeV for 4 years of LHC operation [5J. This result is recognized as the best one 
among others obtained by various methods pQ. 

In the case of a future linear collider (LC) 01 HI E] the most promising method 
for precision M t determination is based on the investigation of the threshold scan of 
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the total cross-section e + e~ — > ti. In this region the form and the height of the cross- 
section are very sensitive to the mass of the top quark. This gives an opportunity 
to determine M t with very high accuracy. A serious difficulty in this approach is a 
precise theoretical calculation of the behavior of the cross-section in the vicinity of the 
threshold, which becomes additionally complicated because of the resonant effects due 
to the strong t—t interaction. A major progress in the calculations was made by way 
of the summation of QCD contributions via solving Lippmann-Schwinger equation 
for the Green function describing the ti production [7j. At the present moment the 
theoretical error of the top mass determined by this method is estimated at 100- 
200 MeV IB1IH], with the experimental error of about 20 MeV [TU]. 

Alternate methods of the M t determination are based on the reconstruction of 
events of the decays of the top quarks. In the basic features they are common at 
LC and at the hadron colliders, but at LC the precision is anticipated better. Thus, 
for example, the systematic error of M t determination by direct reconstruction of ti 
events in e + e~ collisions at y/s = 500 GeV is expected JT] at 340 and 250 MeV in 
hadronic and semi-leptonic channels, respectively, with statistical errors of about 100 
MeV for 1-2 years of LC Since far above the threshold one can expect very 

high precision of the necessary theoretical calculations, the resultant errors should 
be close to that expected at studying the threshold scan of the cross-section. This 
promising anticipation again excites a question about the precision of the top mass 
determination by the method of Ref. jB], but this time in the LC case. Actually this 
method in the LC case has been discussed as a preliminary in (see also review 
0), but the errors have not been determined. So the prospect of this method at LC 
is still not practically known. 

In this article we clear up this question. In contrast to Ref. jB], however, we 
consider the full reconstructed jet of the b quark instead of J/ip or yU + /i~ pair only. 
Such an approach has been considered in [IS], and partially in We follow it by 
keeping in mind that the distribution in any case does emerge in a certain stage 
of the analysis. So from very beginning the analysis can be made in terms of the data 
converted to the form of distribution. (Of course, the systematic errors that arise 
in the course of the converting of the data must be taken into account.) An obvious 
advantage of this approach is a possibility to consider the data in a uniform fashion 
in both cases, LHC and LC. Moreover, this allow us in a simple way to explore a 
possibility of optimization of the algorithm of the extraction of the top mass from 
the data. The elaboration of the latter problem actually is the major purpose of the 
present article. 

In the next section we detail the statement of the problem. In sections |3J and 0] 
we discuss a model for the calculation of the errors. The parameters of the model are 
fixed in section and in the same place the quantitative outcomes are determined. In 
sectionElwe discuss the theoretical uncertainty, and in section[7|we discuss the results. 
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2 Statement of the problem 



We consider the processes 

e + e~ (qq, gg) ^ti^bWbW biv bq x q 2 -> {6-jet + £} + {3 jets} , (1) 

with the 6-jet, isolated lepton i = {e,/i}, and invisible in the final states neutrino 
coming from one t quark, and the remaining three jets coming from another t quark. 
In the experiment the above mentioned states are registered, and measured is a dis- 
tribution 

„, s 1 dcx , . 

Hi) = - -r • 2 

a dq 

Here a is the cross-section of the process (0), q = My is the reconstructed invariant 
mass of the system {6-jet + £}. 

We simulate the results of the experiment under the following suppositions. First 
we suppose that there is a satisfactory method for extracting signal from the data. 
Actually this means the existence of a satisfactory model for the background processes 
that survive after setting of kinematic cuts. 1 Further, we describe the signal in the 
Born approximation, identifying the 6-jet with the b quark. Finally, on the basis of the 
results of Ref. [S] we disregard the effects of finite width of the top quarks. The latter 
assumption means that a _1 da/dq is equal to T^ u dT b £ U /dq, where Y Uv is a partial 
width of the decay t — > biv. (Thus the distribution F becomes process-independent.) 

Direct calculation gives the following formula for the distribution of the partial 
width: 



dq 2 ~ 4V2tt 2 Mf 

arctanf^V-tanf^ 
\M W T W \V W 



(3) 



| (A 2 - q 2 )(q 2 + Ml) + M*,T* W 



Here A 2 = M 2 — Myy, is the total and Tw^iu is the partial width of the W boson, 
and we neglect the masses of the lepton I and the b quark. 2 In this approximation 
IV.^ = 2/9 IV and q ranges between and M t . Fig. [TJshows the distribution F(q) 
defined by formula Q at M f = 170, 175, 180 GeV. From Fig. [T] a dependence of 
F(q) on M t is obvious. So, by comparing the experimental distribution with a set of 
theoretical curves one can determine, in principle, the experimental value of M t . 

In a practical respect, however, it is convenient to compare integrated parameters 
of the distributions. For instance, in Ref. jH] the M t was extracted from the mean 
value (position of the maximum) of the Gaussian distribution approximating the 

1 The set of the cuts and the background processes in the LHC case have been discussed in Ref. jH]- 
In the LC case that has been done in Refs. [Ill 112) . At this stage we do not take manifestly into 
account the kinematic cuts but do that at deriving the quantitative outcomes. 

2 The influence of the mass of the b quark is noticeable at very small g, but this region is inessential 
when considering the moments over the distribution. 
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measured distribution. Refs. f3J E] determined M t by the first moment ( q ) over 
the distribution. In the present article we consider a method of M t determination by 
the higher moments 

pM 

(q n ) = / dq q n F(q) . (4) 
Jo 

Here M is a fixed quantity close to M t , see below for details. In fact this method 
means the matching of the experimental distribution q n F(q) with the corresponding 
theoretical distribution which depends on the parameter M t . 

As we will see below, the insertion of q n factor will significantly increase the 
precision of the M t determination. Eventually this can be checked by a quantitative 
analysis. Nevertheless some hints on this can be a priori seen. Really, with increasing 
n the moment (q n ) becomes in rising measure dependent on the behavior of F(q) in 
a region located between the position of its maximum and a large-g tail where F(q) 
almost vanishes. (More precisely, by the tail we mean a range A < q < M t , where 
in the limit Tw = the distribution identically vanishes by reason of kinematics.) 
Further, in the mentioned region the behavior of F(q) in the greatest measure is 
sensitive to the value of M t , which is seen from Fig. |U As a result, with increasing n 
the sensitivity of (q n ) with respect to M t is increasing. That is why one can expect 
the increasing of the precision of M t extracted from the higher moments. 

Now let us dwell on the details of definition (j3J). The point of the discussion is the 
upper limit in the integral. We set it M instead of conventional M t , meaning the upper 
bound of the region allowed by kinematics, because M t also is a parameter which is 
subject to determination. In order to avoid an inconvenience, we put in the place of 
the upper limit a certain predetermined value M fixed close to M t . Simultaneously 
we adjust the normalization of the distribution F(q) so that to provide the equality 
(1) = 1. The moments (q n ) at n > 1 after this redefinition practically do not change 
(at not too large n) in view of almost vanishing F in the tail at large q. 

So, we define the experimentally measured value of the top quark mass as a 
solution to the equation 

(?"> = (<f>exp • (5) 

Here in the r.h.s. the moment is determined (at a given M) on the basis of the 
experimental data, and that in the l.h.s. on the basis of the theoretical distribution 
which depends on the parameter M t . Let at a given n a solution to the equation © 
be M t = M t ( n \. Then the error of the solution can be determined as 

(6) 

M t =M t(n) 

Our aim is to estimate AM t ( n ) and find an optimal value of n which would minimize 
AM t ( n ). Since by virtue of © the derivative d(q n )/dM t is known, the problem is 
reduced to the determination of the statistical and systematic errors, the components 
of the experimental error A(g n ) exp . 



AM, 



t(n) 



A<0 



exp 



d(q n ) 
diWT 
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3 Statistical errors 



We determine the statistical errors of the moments on the supposition that the data 
averaged over ensemble are described by F(q) = T^dT^/dq with T Uv determined 
by formula © at M t = 175 GeV. 

Let 5qi be the size of a bin, within which z-th element of the distribution is 
measured, and let iVj be the number of events counted in this bin on the average. 
Then 

F(q i )6q i = W i /N. (7) 

Here N is the total number of events counted in all bins on the average. Further we 
do not distinguish between N and N = ^ Ni, the total number of events counted in 
all bins in the given experiment. The experimentally measured n-th moment is 

N- 

(<f>e XP = E^- (8) 

i 

By virtue of (JJJ) the averaged experimental moment (q n ) exp is found by formula (j3J). 
Since iVj is distributed by Poisson law with parameter Ni, the variance of {q n ) e xp is 

i 

Formula © implies the following estimation for the statistical error: 



A sta V) exp = 




To give an idea of the behavior of A stat (g n ) exp , we present in Fig.Elby a continuous 
curve the ratio A stat (g ri ) exp / (<? n )ex P calculated at N = 4000 (corresponds to LHC case, 
see Sect. EJ. It is seen from the figure that with increasing n the ratio is growing. 
This is explained by the shift (to the right) of the position of maximum of q n F(q) 
from the position of maximum of F(q), where the statistics is largest. As a result a 
statistical reliability of (q n ) e xp comes down. Another important property of the ratio 
is the change of the mode of the growth beginning with n ~ 15. This is explained by 
the emergence of a noticeable contribution from the large-g tail in q n F(q). The latter 
property is illustrated by the set of the curves represented by Fig. El 

In fact the emergence of a noticeable contribution from the tail is an undesirable 
effect since on the tail the uncertainty from the background is comparable with the 
signal process. In order to avoid this difficulty one can correct the definition of the 
moments by introducing a cutoff in the integral in (J3J). The position of the cutoff 
should be determined so that to isolate the second (unphysical) peak in the tail 
of q n F(q) but simultaneously to keep as much as possible a statistical significance 
of the sample of events. It is clear that the optimal cutoff should be placed in the 
neighborhood of a local minimum between the two peaks of q n F(q) (if the second peak 
appears). From Fig. [Slit is seen that at n ~ 40 the sought-for point is distant by about 
two half- widths to the right of the position of the maximum of q n F(q) . So a simplified 
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algorithm for the cutoff may be determined by setting A n = min{g nextr + 2r rir i g ht, M}, 
where g nex tr is the position of the maximum of q n F(q) and r nright is the half-width 
from the right. Thus we come to the following definition of the effective moments: 

(q n V S = j\ Q n F(q) j f\ F{q) . (11) 

In the experimentally determined effective moments the cutoff must be the same. 
Ultimately A stat (g")^ p is defined by formula (fTUj) with (q 2n ) replaced by (g 2n ) eff but 
with introducing the cutoff A n instead of A2 n - The latter anomalous prescription 
follows immediately from the derivation of formula (jlUj) . 

The behavior of A stat (g n )^ p /(g")^L is shown by the dashed curve in Fig. El It 
is seen from the figure that the transition to the effective moments implies no no- 
ticeable modification up to n « 15, while at the larger n the growth of the ratio 
becomes stabilized. A similar behavior is observed in the basic formalism (without 
the transition to the effective moments) in the limit Ty/ — > 0, when the large-g tail 
identically vanishes. 



4 Systematic errors 

Proceeding to the systematic errors it is necessary at first to clarify a reason of 
their origin. For this purpose we use the analysis of Ref. jH] of the errors of the 
MjM distribution simulated with the PYTHIA and/or HERVIG event generators. 
By the main sources of the systematic errors Ref. |)J found the uncertainties in the 
b quark fragmentation (including the final state radiation) and the uncertainties in 
the background processes. It is clear that the same sources should be the main ones 
at solving the inverse problem, the determination of the Mu distribution from Mjjyt 
distribution which is considered virtually as the data. In the LC case we expect the 
same pattern of the origin of systematic errors. 

On this basis we consider at first the error resulting from the uncertainty in the b 
quark fragmentation. For brevity we call it by the error of the type I. At the level of 
Mb£ distribution it appears as the uncertainty in the number of the bin within which 
the number of events, Ni, is measured. In the continuous case this error becomes the 
uncertainty Aq in the determining of q variable. 

Suppose that Aq is sufficiently small. Then, neglecting the nonlinear effects, we 
have 

A^ s V)exp= / dg [q n F(q)]' Aq. (12) 
Jo 

Here the prime means the derivative with respect to q. The systematic error I of the 
effective moment (q n )t^p ^ s estimated similarly, with replacing the upper bound M 
by A n and, then, dividing the result by the normalization factor as in formula (jllj) . 
The normalization factor itself should be the same as it controls the total number of 
events that are taken into consideration at the determining of the effective moment. 

The determination of Aq we carry out with the aid of the following reasoning. 
First we note that the invariant mass q 2 actually is the doubled scalar product of 



6 



4-momenta of the b quark and the lepton I. So in the laboratory frame it can be 
represented as q 2 = E^K, where is the energy of the b quark, and if is a factor 
proportional to the energy of the lepton £. (Additionally K includes a dependence 
on angular variables which, however, is relatively weak.) Further, by the calculating 
of the differential we get Aq = | (AEb/Eb + AK/K) q, where AEb and AK are the 
corresponding errors. A more precise estimation is determined by the sum in the 
quadratures. Thus we come to a linear dependence with a certain coefficient r, 



Aq = r q , r = — 



AE b \ 2 . fAK\ 2 



1 1/2 

(13) 



The systematic error arising after subtraction of the background we name by the 
error of type II. It appears in the absolute value of the distribution function. So it 
should be described as an additive contribution SF to function F. Correspondingly, 
we get the following formula for the error II of the moments: 

fAI 

A sysI V) cxp = / dg q n 5F(q). (14) 
Jo 

The error II of the effective moments {q n )l^ p is defined by a similar formula to within 
modifications listed below 

It is reasonable to determine SF(q) on supposition that it vanishes at the bound- 
aries of the phase space and when passing from small q to large q it only once changes 
the sign. The simplest form of a function satisfying to these requirements is a poly- 
nomial of degree three, 

5F = hq(q — M/2) (q — M). (15) 

Parameter h in ()15|) describes the amplitude of the error and it is subject to further 
determination. 



5 Numerical results 

We assign the following values for the parameters having a global meaning: 

M w = 80.4 GeV, IV = 2.1 GeV, M t = M= 175 GeV. (16) 

The remaining parameters N, r, and h depend on the conditions of the consideration. 
Recall that N means the volume of the representative sample of events, parameter 
r characterizes the error in the invariant mass of the M-sjstem, and h describes the 
error arising after the subtraction of the background processes. 

With reference to LHC case, the parameters N, r, h we fixe on the basis of 
the results of Ref. j^j. Since in that work the Mj^pe distribution was determined at 
iV = 4000 (with kinematic cuts and for 4 years of LHC), in our investigation we 
set this value for N, as well. Parameters r and h we fix based on the properties 
of Mjfyi distribution and the direct results derived in Ref. jB] from these properties. 
First we use the estimation A sys (Mj/^) = +0.3/— 0.4 GeV and the derived from it 
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n 


A stat M t(n) 


A sys l M t{n) 


A sys n M t{n) 


A sys M t(n) 


AM f(n) 


1 


2.07 


0.62 


0.30 


0.69 


2.18 


5 


0.62 


0.13 


0.20 (0.17) 


0.24 (0.21) 


0.66 


10 


0.45 


0.07 (0.06) 


0.20 (0.13) 


0.21 (0.14) 


0.50 (0.48) 


15 


0.41 (0.40) 


0.05 (0.03) 


0.24 (0.12) 


0.24 (0.12) 


0.48 (0.42) 


20 


0.42 (0.39) 


0.03 (0.00) 


0.32 (0.11) 


0.32 (0.11) 


0.52 (0.40) 


30 


0.59 (0.38) 


0.02 (0.03) 


0.63 (0.11) 


0.63 (0.11) 


0.86 (0.39) 



Table 1: Statistical, systematic I and II, and the systematic summed in quadrature errors 
presented in GeV in the LHC case. The last column represents the sum of the statistical 
and the systematic errors. In the brackets we show the results calculated by the method of 
the effective moments (if they are different with the results calculated by the basic method) . 

the result AM t = +0.6/— 0.8 GeV. Considering in the framework of our investigation 
the latter quantity as the uncertainty of the input parameter M t , we get by direct 
calculation A sys (g) exp = +0.47/- 0.62 GeV. By comparing this with A sys (M m ) we 
obtain an energy scale factor of 1.6, which describes the spreading of the M^g distri- 
bution when converting it to the Mm distribution. Using the mentioned factor, from 
A sysn (M JAW ) exp < 0.15 GeV © we further derive an estimation A sysII (g) exp < 0.24 
GeV. From this result and (fT3j) we get h ~ 1.7 x 10 -10 GeV -4 . (Hereinafter we 
take the upper bounds as the estimations.) 

Knowing A sys (g) cxp and A sys n (g) exp , we immediately get A sysI (g) exp ~ +0.41/— 
0.57. From here and formula (fT2*|) there follows r ~ 0.004-0.006. Further we use 
the average value r = 0.005. It is worth noticing that the same estimation for r 
follows from formula (|T3"j) when taking into consideration the 1%-precision of the 
determination of the energy of b jets expected at LHC [T] , and additionally neglecting 
AK/K as compared to AE b /E b . 

Now as we know N, r, h, we can calculate A stat (g n ) exp and A sys I,n (g n ) e xp at any n. 
Then we calculate A stat M^ n ) and A sys I,n M t ( n ) . The dependence on n of these errors 
is shown by the solid lines on Figs. EHH1 The dashed lines show the errors obtained 
by the method of the effective moments. (The break of slope in the dashed line in 
Fig. 5 is explained by the change of the sign in the integral in formula ()12j) appearing 
after introducing the cutoff A n .) In Tabled we present the numerical results at some 
n. In the same place we show the summed in the quadrature errors A sys M t ( n ) and 
AM t ( n ). It should be noted that at n = 1 the systematic errors in Tabled practically 
coincide with those in jHj. The reason is that we have fixed the parameters of the 
model actually by matching the errors of the first moments. 

In the LC case, unfortunately, there are no published results that could allow us 
in a similar way to fix the parameters of the model. Therefore we make use mainly 
of indirect methods. Parameter iV we fix by the following reasoning. First we note 
that er(e + e~ — > tt) ~ 0.6 pb at s/s = 500 GeV. So at the integrated luminosity of 300 
fb -1 , corresponding to 1-2 years of running, approximately 180 000 tt pairs must be 
generated. Since the branching of the process (0) is near of 30%, only 54 000 events 
of tt are related to our investigation. The efficiency of their detection we estimate as 
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n 


A stat M t(n) 


A sys l M t{n) 


A sys n M t(n) 


A sys M i(n) 


AM f(n) 


1 


0.80 


0.17 


0.19 


0.26 


0.84 


5 


0.24 


0.04 


0.13 (0.11) 


0.13 (0.12) 


0.27 


10 


0.17 


0.02 


0.13 (0.09) 


0.13 (0.09) 


0.27 (0.20) 


15 


0.16 


0.01 


0.15 (0.08) 


0.15 (0.08) 


0.22 (0.17) 


20 


0.16 (0.15) 


0.01 (0.00) 


0.21 (0.07) 


0.21 (0.07) 


0.26 (0.16) 


30 


0.23 (0.15) 


0.01 


0.41 (0.07) 


0.41 (0.07) 


0.46 (0.16) 




Table 2 


: The same 


that in Tabled 


in the LC case. 



follows. Suppose that at LC the efficiency of the detecting of W^-pairs decaying in 
a semi-leptonic channel will be the same as at LEP2, i.e. ~ 80% [T^j. In addition, 
following Ref. we suppose that the 6-jet tagging efficiency at LC will be about 
80%. In summary this gives an acceptance of 50% which implies N = 27000. 

Parameter r we fix based on the systematic error of M t obtained in ^T] in the 
approach of the direct reconstruction of tt events in the semi-leptonic channel. Addi- 
tionally we use the note that in the kinematic range near the upper Mm endpoint the 
neutrino practically does not contribute to the total invariant mass of the decay prod- 
ucts of the top quark. Therefore the determination of the M t in the mentioned range 
practically is the same that the determination of the Mm invariant mass. Ref. fl] 
obtained A sys M t = 250 MeV. So we set A sys M M = 250 MeV. Assuming that this is 
the error of the type I, we equate Aq to this value. Finally, by setting Aq = rq, 
q ~ M t we get r = 0.0014. 

Parameter h we fix by proceeding to the note of Ref. [Oj about the decreasing 
of the systematic error II of the average (Mjm) below of 0.1 GeV at the increasing 
of statistics up to iV ~ 10 4 . From this by using the above method we get a rough 
estimate h~l.lx 10~ 10 GeV" 4 . 

Knowing N, r, h, we find A stat (g n ) exp , A sys I ' II (g n ) exp and then A stat M t(n ) and 
A sys l,ll Mt( n ). Since in our model the pattern of the dependence on n is the common 
one in the LC and LHC cases, the difference between these cases appears in the scales 
of the errors only. This allows us to present the results on Figs. EHSlwith adding new 
scales. The numerical results are presented in Table 121 It is interesting to note that 
A sys I M t ( 1 ) turns out smaller than A sys M t obtained in Ref. ^1] in the framework of 
the direct reconstruction of events. Nevertheless this does not mean an inconsistency. 
Really, we equate A sys M t of ^T] to Ag at q ~ M t , but the dominant contributions to 
A sys I (g) CX p are formed at strictly smaller q than M t , which is obvious from formulas 
(JHJ) and fEIJ). This effect diminishes A sys I M t{1) compared to A sys M t of [Hj. 

6 Theoretical uncertainty 

The analysis of the previous sections shows that the experimental accuracy of the 
M t determination can be considerably improved by the transiting to the high degrees 
of the moments. So, for the achieving of the eventual high accuracy a theoretical 



9 



uncertainty becomes more and more crucial. In this connection it is important to 
understand whether the theoretical error in the M t determination can be made smaller 
than the experimental error at the high degrees of the moments. If it will be possible 
then the theoretical error will not spoil the expected accuracy. Below we discuss this 
question in somewhat qualitative manner since the possibility of the solution first of 
all is important. 

Let us begin with the note that the origin of the theoretical uncertainty in the M t 
determination is connected with the uncertainty of the calculation of the theoretical 
moment (q n ) in equation (JHJ). Further, at the determining of AM t t n \ the corresponding 
error A th {q n ) is to be added (in quadratures) to A(q n ) exp in formula ©. So the 
problem is reduced to the question about a possibility of carrying out the calculations 
so precisely that to keep the error A th (g n ) smaller than A(g n ) eX p. 

In practice it is convenient to compare relative errors like A{q n ) / (q n ) instead of 
the proper errors A(q n ). Fortunately the experimental relative error A(q n ) cxp / (q n ) e xp 
is growing at the transiting to the high degrees of the moments; its behavior is similar 
to that represented in Fig. El In particular, at transiting from n = 1 to n — 15 the 
/ ^(? n )o xp /(q ,n )exp increases from 1.8% to 5.2%(4.5%) in the LHC case, and from 0.7% 
to 2.4%(1.9%) in the LC case. So with increasing n the requirement for the theoretical 
relative error A th (q n ) / (q n ) is weakening. 

Generally a theoretical error arises from a parametric uncertainty and an intrinsic 
uncertainty of the itself calculation. The parametric uncertainty originates mainly 
from the parameters that are worse known. In the given case they are the widths of 
the W boson and of the top quark. The analysis of Ref . |H] shows that the uncertainties 
in these parameters practically do not affect the first moment. Moreover, even the 
switching-off of the widths gives a negligible effect. In the 15-th moment, the varying 
of Yw within the experimental error ATw = 0.04 GeV results in A(g 15 )/(g 15 ) = 
0.09%(0.06%), which is insignificant, as well. Unfortunately we cannot estimate the 
variance of the moments with varying the width of the top quark, since from the very 
beginning we use for the top quarks the narrow width approximation. Nevertheless 
by basing on the results of Ref. [U] we expect a negligible variance of the moments 
in this case too. This is corroborated by the lack of reasons leading to appreciably 
greater sensitivity of the moments with respect to the width of the top quarks than 
to the width of the W boson. 

It should be mentioned, however, that the complete switching-off of the widths 
can vary noticeably the high-degree moments. Thus, the setting Tyy = implies a 
shift of (g 15 ) on 4.6%(3.2%) which can be compared with the above estimations for 
the experimental relative errors. This means that the calculation of the high-degree 
moments must be carried out with the taking into consideration of the realistic values 
of the widths. The latter requirement, of course, is unnecessary for the estimation of 
the errors only, the case of the investigation of the present article. 

Now let us consider the errors of the itself calculation. First we note that all the 
processes in ((TJ) go far above the thresholds of the production of unstable particles, the 
W bosons and the top quarks. Therefore their production and decay can be described 
by the standard methods [TEj . namely with the Dyson resummation in the leading 
order calculation and in the pole approximation at calculating the perturbation-theory 
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corrections. 3 Thus, the problem is reduced to the estimation of the order of the 
perturbation theory, which is necessary for satisfying the required precision of the 
calculation. 

Further we note that the corrections to the moments are to be calculated by way 
of calculating the corrections to the distribution F(q). For kinematic reasons the 
behavior of the latter corrections in the basic features should follow the behavior of 
the distribution. Namely, A th F(q) must vanish at the ends of the kinematic region 
because of the vanishing phase volume. Furthermore, A th F(q) must almost vanish 
on the tail at large q due to the smallness of the width of the W boson. (Recall that 
in the limit = the distribution is completely suppressed on the tail by reason 
of kinematics.) So A th F(q) must be precisely known mainly in the middle of the 
kinematic region but not near its ends including the tail. At transiting to the high- 
degree moments this condition is maintained. Moreover, in some sense it becomes 
even more strong. Really, at low q the contributions to the moments are additionally 
suppressed by the factor q n . At large q, in the case of the effective moments, the 
contributions are completely suppressed by the cutoff A n . In addition, the larger 
n the more the distant between the cutoff and the M t , the right boundary of the 
actual range of kinematic variable (since A n — > A from the right as n — ► oo). In 
particular, Ai = 171 GeV but A 15 = 160 GeV which is distant form M t by 15 GeV. 
The mentioned feature is valuable for our consideration as the cut-off of the ends of 
the region of the kinematic variable implies a suppression of large logarithms that can 
arise near the ends at calculating the perturbation-theory corrections. In the final 
analysis this allows us to use the naive counting method for estimating the corrections 
to the moments. 4 

With this in mind we farther use a rather rough approach which is based on a 
comparison between the corrections to the moments and to the width of the top quark. 
(Notice that the width actually is the zero moment accurate to the normalization.) 
The key reason of the approach is the observation that the integrals for the moments 
and the width, and for the corrections to the moments and the width, accumulate the 
contributions mainly from the middle region of the kinematic variable. So, supposing 
that in this region the correction to the distribution A th F(q) varies weakly in the 
units of F(q), one can expect that the corrections to the moments and to the width 
should be close to each other in the relative units. By the closeness we admit here 
a factor of order of several units. It is worth mentioning that even in the case of 
the experimental errors, which depend strongly on the shape of the distribution, the 
relative errors at n — 1 and n = 15 differ from each other by a factor of 2.5-3.5 only. 

As we know, the electroweak one-loop correction to the top quark width amounts 
approximately 2%. The QCD one-loop correction is near of 10%, while the two-loop 
one is near of 2%. (See PQ and the references therein.) The comparison of these values 

3 As variants, one can exploit the method of an effective field theory for calculating the resonant 
processes ^7] or the modified perturbation theory based on distribution theory [THIH^] . 

4 It should be emphasized that we discuss here the corrections to the t — > b transition but not to 
the 6-quark fragmentation, including the perturbative fragmentation. The latter process is described 
by the convolution of the cross-section with the fragmentation function, and this operation is to be 
fulfilled in the framework of MC event generators. 
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with the above estimations for the experimental relative errors demonstrates that the 
one-loop electroweak and two-loop QCD corrections are enough to remain within the 
required limits (in both cases, LHC and LC). It should be noted that the mentioned 
corrections to the distribution can be certainly calculated since the corrections to the 
width have been calculated. Finally we note also that only the direct calculations 
of these corrections can explicitly solve the problem of the theoretical errors to the 
moments. 

The mentioned calculation, however, would not yet entirely close the problem of 
the theoretical error of the M t determination because of the problem of nonpertur- 
bative nature caused by a renormalon contribution. Below for the completeness we 
only briefly consider this problem as its solution is known, at least in a conceptual 
respect. The problem actually is connected with the kind of the mass which is to be 
determined through an experimental measurement. In fact there are different masses, 
but only a Lagrangian mass is ultimately valuable since only the Lagrangian mass 
can be constrained with other fundamental parameters of the theory. The important 
representatives of the Lagrangian mass are the pole and MS masses. The directly 
measurable one is the pole mass, which is determined by kinematics. Correspond- 
ingly, the currently used algorithms of extracting M t from the data are turned to 
the pole mass. However, because of the renormalon contribution the pole mass de- 
termination faces an extra uncertainty of order of O(Aqcd) [20J. Numerically it can 
amount hundreds of Me Vs. 

The above mentioned difficulty can be bypassed in the framework of the following 
algorithm (below we state one of its possible variants) j20j. First, all theoretical 
calculations are to be fulfilled in the terms of the pole mass. Then the value of the 
pole mass is to be determined from the matching with data. Remember, at this stage 
the result includes the renormalon contribution. Further, by means of the well-known 
formula relating the pole mass with the MS mass (see Ref. for example), the MS 
mass is determined. At this step the result gets again the renormalon contribution 
but, as is declared, it cancels the previous one. (So the inaccuracy in the relation 
between the pole and MS masses is charged to the pole mass.) Direct calculations in 
certain examples [HI E] demonstrate effectiveness of the above algorithm. 

So, the problem is initially stated as though for the pole mass determination, but 
at the final stage the MS mass is determined. This allows one to avoid a theoretical 
systematic uncertainty of order of O(Aqcd) caused by the renormalon contribution. 
Returning to our outcomes, we see that the theoretical error of the top mass deter- 
mination can be really made smaller than the experimental error. 

7 Discussion 

The major result of this article is the detection of the effect of decreasing of statistical 
and systematic errors of the top quark mass measured from Mu distribution, when 
applying the technique of the moments and proceeding to the high degrees of the 
moments. The optimal value of the degree minimizing the errors is found near n = 15. 
For the determining of the errors we have attracted a simple enough model. Its 



12 



parameters in the LHC case have been fixed on the basis of the results obtained 
earlier [H] by the MC modelling method. As applied to LC the parameters have been 
fixed mainly by indirect methods. Knowing the parameters and the dependence on 
the degree n of the moments, we have estimated the errors at varied n and have 
found the optimal value of n, minimizing the errors. The optimal value n = 15 is 
clearly visible in the framework of the basic method of calculating the moments. The 
applying of the technique of the effective moments diminishes the errors at n — 15 
by 10-20%, but at the further increasing of n the results practically do not vary (see 
Figs. HE1 and Tables EH2J) • 

At the optimal value n = 15 the total error AM t is found close to 500 MeV in the 
LHC case, and close to 200 MeV in the LC case. In the LHC case the above accuracy 
more than in twice exceeds the accuracy obtained by the other methods pQ including 
the original method of Ref. jB]. In the LC case the estimated accuracy of the M t 
determination is close to that expected at scanning the ti production threshold jHUHl- 

In conclusion it should be mentioned, once again, that at the intermediate stage 
of the analysis we have introduced simplifications allowing to minimize calculations. 
However at the final stage all estimations have been made on the basis of realistic 
values of the parameters. This peculiarity should not reduce the legitimacy of the 
detected behavior of the errors and, moreover, of their rounded estimations. Never- 
theless the quantitative outcomes could be improved by further calculations based on 
the direct applying of a proper MC event generator. 
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Figure 1: Distribution F(q) = T^dT/dq, 
q = Mm, at M t = 170, 175, and 180 GeV. 
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Figure 2: The ratio A stat (^) cxp /(g m ) cxp de- 
pending on n {M t = 175 GeV, N = 4000). 
The continuous curve represents the results 
described by formula (|10jl. The dashed 
curve represents the results obtained by the 
method of the effective moments. 



Figure 4: The statistical error A stat M t(n ) 
depending on n. The dashed curve repre- 
sents the results obtained by the method of 
the effective moments. The left and right 
vertical axes scale in GeV the results for LHC 
and LC cases, respectively. 
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Figure 5: The same that on Fig. 21 for 

A sysI M, (n) . 
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Figure 3: The shape of the function q n F(q) 
at M t = 175 GeV, n = 1, 5, 15, 40 (in arbi- 
trary normalization). 
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Figure 6: The same that on Fig. |1] for 

a^ s n M t(n) . 
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